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!>■ ! 1 Introduction 

Radio regulatory bodies are recently recognizing that rigid spectrum assignment granting exclusive use to 
c/3 , licensed services is highly inefficient., due to the high variability of the traffic statistics across time, space, 
l— ~~ '> and frequency. Recent Federal Communications Commission (FCC) measurements show that, in fact, 
the spectrum usage is typically concentrated over certain portions of the spectrum, while a significant 

^ ! 

amount of the licensed bands (or idle slots in static time division multiple access systems with bursty 

C — . 

Q^ ■ traffic) remains unused or underutilized for ninety percent of time [1]. It is not surprising then that 
O " 

this inefficiency is motivating a flurry of research activities in engineering, economics and regulation 

' communities in the effort of finding more efficient spectrum management policies. 
OO ! 

As pointed out in many recent works [2l EI [H [5] , the most appropriate approach to tackle the great 
spectrum variability as a function of time and space calls for dynamic access strategies that adapt to 
the electromagnetic environment. Cognitive Radio (CR) originated as a possible solution to this problem 
[6] obtained by endowing the radio nodes with "cognitive capabilities", e.g., the ability to sense the 
electromagnetic environment, make short term predictions, and react intelligently in order to optimize 
the usage of the available resources. Multiple paradigms associated with CR have been proposed (SUHJHIE], 
depending on the policy to be followed with respect to the licensed users, i.e. the users who have acquired 
the right to transmit over specific portions of the spectrum buying the relative license. The most common 
strategies adopt a hierarchical access structure, distinguishing between primary users, or legacy spectrum 
holders, and secondary users, who access the licensed spectrum dynamically, under the constraint of not 
inducing Quality of Service (QoS) degradations intolerable to the primary users. Within this context, 
three basic approaches have been considered to allow concurrent communications: spectrum overlay, 



underlay and interweave 



1 There is no strict consensus on some of the basic terminology in cognitive systems [4]. Here we use interweave as in [5] 
which is sometimes referred to as overlay communications [4j. 
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In overlay systems, as proposed in [7], secondary users allocate part of their power for secondary 
transmissions and the remainder to assist (relay) the primary transmissions. By exploiting sophisticated 
coding techniques such as dirty paper coding based on the knowledge of the primary users' message and/or 
codebook at the cognitive transmitter, these systems offer the possibility of concurrent transmissions 
without capacity penalties. However, although interesting from an information theoretic perspective, 
these techniques are difficult to implement as they require noncausal knowledge of the primary signals at 
the cognitive transmitters. 

In underlay systems, the secondary users are also allowed to share resources with the primary users, 
but without any knowledge about the primary users' signals and under the strict constraint that the 
spectral density of their transmitted signals fall below the noise floor at the primary receivers. This 
interference constraint can be met using spread spectrum or ultra-wideband communications from the 
secondary users. Both transmission techniques do not require the estimation of the electromagnetic 
environment from secondary users, but they are mostly appropriate for short distance communications, 
because of the strong constraints imposed on the maximum power radiated by the secondary users. 

Conversely, interweave communications, initially envisioned in [6], are based on an opportunistic or 
adaptive usage of the spectrum, as a function of its real utilization. Secondary users are allowed to adapt 
their power allocation as a function of time and frequency, depending on what they are able to sense 
and learn from the environment, in a nonintrusive manner. Rather than imposing a severe constraint on 
their transmit power spectral density, in interweave systems, the secondary users have to figure out when 
and where to transmit. Differently from underlay systems, this opportunistic spectrum access requires an 
opportunity identification phase, through spectrum sensing, followed by an opportunity exploitation mode 
[4]. For a fascinating motivation and discussion of the signal processing challenges faced in interweave 
cognitive radio systems, we suggest the interested reader to refer to [2]. 

In this paper we focus on opportunistic resource allocation techniques in hierarchical cognitive net- 
works, as they seem to be the most suitable for the current spectrum management policies and legacy 
wireless systems [4]. We are specifically interested in devising the most appropriate form of concurrent 
communications of cognitive users competing over the physical resources let available from primary users. 
Looking at opportunistic communication paradigm from a broad signal processing perspective, the sec- 
ondary users are allowed to transmit over a multi-dimensional space, whose coordinates represent time 
slots, frequency bins and (possibly) angles, and their goal is to find out the most appropriate transmission 
strategy, assuming a given power budget at each node, exploring all available degrees of freedom, under 
the constraint of inducing a limited interference, or no interference at all, at the primary users. 

In general, the optimization of the transmission strategies requires the presence of a central node 
having full knowledge of all the channels and interference structure at every receiver. But this poses a 
serious implementation problem in terms of scalability and amount of signaling to be exchanged among 
the nodes. The required extra signaling could, in the end, jeopardize the promise for higher efficiency. 
To overcome this difficulty, we concentrate on decentralized strategies, where the cognitive users are able 
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to self-enforce the negotiated agreements on the spectrum usage without the intervention of a centralized 
authority. The philosophy underlying this approach is a competitive optimality criterion, as every user 
aims for the transmission strategy that unilaterally maximizes his own payoff function. The presence of 
concurrent secondary users competing over the same resources adds dynamics to the system, as every 
secondary user will dynamically react to the strategies adopted by the other users. The main question 
is then to establish whether, and under what conditions, the overall system can eventually converge 
to an equilibrium from which every user is not willing to unilaterally move, as this would determine a 
performance loss. This form of equilibrium coincides with the well-known concept of Nash Equilibrium 
(NE) in game theory (see, e.g., [HE]). In fact, game theory is the natural tool to devise decentralized 
strategies allowing the secondary users to find out their best response to any given channel and interference 
scenario and to derive the conditions for the existence and uniqueness of NE. 

Within this context, in this paper, we propose and analyze a totally decentralized approach to design 
cognitive MIMO transceivers, satisfying a competitive optimality criterion, based on the achievement 
of Nash equilibria. To take full advantage of all the opportunities offered by wireless communications, 
we assume a fairly general MIMO setup, where the multiple channels may be frequency channels (as 
in OFDM systems) [H}]-[T2], time slots (as in TDMA systems) pUJdl], and/or spatial channels (as in 
transmit/receive beamforming systems) [H]. Whenever available, multiple antennas at the secondary 
transmitters could be used, for example, to put nulls in the antenna radiation pattern of secondary 
transmitters along the directions identifying the primary receivers, thus enabling the share of frequency 
and time resources with no additional interference. Our initial goal is to provide conditions for the 
existence and uniqueness of NE points in a game where secondary users compete against each other to 
maximize their performance, under the constraint on the maximum (or null) interference induced on the 
primary users. The next step is then to describe low-complex totally distributed techniques able to reach 
the equilibrium points of the proposed games, with no coordination among the secondary users. 

2 System Model: Cognitive Radio Networks 

We consider a scenario composed by heterogeneous wireless systems (primary and secondary users), as 
illustrated in Figured! The setup may include peer-to-peer links, multiple access, or broadcast channels. 
The systems coexisting in the network do not have a common goal and do not cooperate with each other. 
Moreover, no centralized authority is assumed to handle the network access from secondary users. Thus, 
the secondary users are allowed, in principle, to compete for the same physical resources, e.g., time, 
frequency, and space. We are interested in finding the optimal transmission strategy for the secondary 
users, using a decentralized approach. A fairly general system model to describe the signals received by 
the secondary users is the Gaussian vector interference channel: 




(1) 
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Figure 1: Hierarchical cognitive radio network with primary and secondary users. 

where x q is the riT q -dimensional block of data transmitted by source q, H qq is the riR q x riT q (complex) 
channel matrix between the q-th transmitter and its intended receiver, H rg is the hr x ny r cross-channel 
matrix between source r and destination q, y q is the nR q - dimensional vector received by destination q, and 
n q is the n# ? -dimensional noise plus interference vector. The first term in the right-hand side of (pQ) is the 
useful signal for link q, the second and third terms represent the Multi-User Interference (MUI) received 
by secondary user q and caused from the other secondary users and the primary users, respectively. The 
vector n q is assumed to be zero-mean circularly symmetric complex Gaussian with arbitrary (nonsingular) 
covariance matrix R„ 9 . For the sake of simplicity and lack of space, we consider here only the case where 
the channel matrices H qq are square nonsingular. We assume that each receiver is able to estimate the 
channel from its intended transmitter and the overall MUI covariance matrix (alternatively, to make 
short term predictions, with negligible error) The receiver sends then this information back to the 
transmitter through a low bit rate (error-free) feedback channel, to allow the transmitter to compute the 
optimal transmission strategy over its own link. 

The model in ([T]) represents a fairly general MIMO setup, describing multiuser transmissions over 
multiple channels, which may represent frequency channels (as in OFDM systems) [H)]-[12], time slots (as 
in TDMA systems) pUKHJ, or spatial channels (as in transmit/receive beamforming systems) [13J. Dif- 
ferently from traditional static or centralized spectrum assignment, the cognitive radio paradigm enables 
secondary users to transmit with overlapping spectrum and/or coverage with primary users, provided 
that the degradation induced on the primary users' performance is null or tolerable. How to impose inter- 
ference constraints on secondary users is a complex and open regulatory issue [21 0]. Roughly speaking, 
restrictive constraints may marginalize the potential gains offered by the dynamic resource assignment 
mechanism, whereas loose constraints may affect the compatibility with legacy systems. Both determinis- 
tic and probabilistic interference constraints have been suggested in the literature [H El SI HI] , namely: the 

2 How to obtain both channel-state information and MUI covariance matrix estimation goes beyond the scope of this 
paper; the interested reader may refer to, e.g., PUS], where classical signal processing estimation techniques are properly 
modified to be successfully applied in a cognitive radio environment. 
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maximum MUI interference power level perceived by any active primary user (the so-called interference 
temperature limit) [TJ [2] and the maximum probability that the MUI interference level at each primary 
user's receiver may exceed a prescribed threshold [H,[T5]. In the presence of sensing errors, the access to 
channels identified as idle should also depend on the goodness of the channel estimation. As shown in 
[T7] . in this case the optimal strategy is probabilistic, with an probability depending on both the false 
alarm and miss probabilities. 

In this paper we are primarily interested in analyzing the contention among the secondary users over 
a multiuser channel where there are primary users as well. To limit the complexity of the problem, in the 
effort to find out distributed techniques guaranteed to converge to NE points, we restrict our analysis to 
consider only deterministic interference constraints, albeit expressed in a very general form. In particular, 
we envisage the use of the following possible interference constraints (see also Figure [2]): 

Co.l Maximum transmit power for each transmitter: 

s{\\x q \\l}=Tr(Q q )<P q , (2) 

where Q q denotes the covariance matrix of the symbols transmitted by user q and P q is the transmit 
power in units of energy per transmission. 

Co. 2 Null constraints: 

Uf Q q = 0, (3) 

where XJ q is a strict tall matrix (to avoid the trivial solution Q q = 0), whose columns represent the spatial 
and/or the frequency "directions" along with user q is not allowed to transmit. We assume, without loss 
of generality (w.l.o.g.), that each matrix XJ q is full-column rank. 

Co. 3 Soft shaping constraints: 

Tr (G q Q q G q ) < P q ve , (4) 

where the matrices G q are such that their range space identifies the subspace where the interference level 
should be kept under the required thresholds 

Co. 4 Peak power constraints: the average peak power of each user q can be controlled by constraining the 
maximum eigenvalue [denoted by A max (-)] of the transmit covariance matrix along the directions spanned 
by the column space of G q : 

A max (Gf Q.G,) < pP eak , (5) 

where pP eak {$ the maximum peak power that can be transmitted along the spatial and/or the frequency 
directions spanned by the column space of G q . 



3 The interference temperature limit constraint [2] is given by the aggregated interference induced by all secondary users. 
In this paper, we assume that the primary user imposing the soft constraint, has already computed the maximum tolerable 
interference power Pg VC for each secondary user. The power limit Pq Ve can also be the result of a negotiation or opportunistic 
based procedure between primary users (or regulatory agencies) and secondary users. 
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Figure 2: Example of null/soft shaping constraints. 

The structure of the null constraints in j3]) is a very general form to express the strict limitation 
imposed on secondary users to prevent them from transmitting over the subchannels occupied by the 
primary users. These subchannels are modeled as vectors belonging to the subspace spanned by the 
columns of each matrix XJ g . This form includes, as particular cases, the imposition of nulls over: 1) the 
frequency bands occupied by the primary receivers; 2) the time slots occupied by the primary users; 3) 
the angular directions identifying the primary receivers as observed from the secondary transmitters. In 
the first case, the subspace is spanned by a set of IFFT vectors, in the second case by a set of canonical 
vectors, and in the third case by the set of steering vectors representing the directions of the primary 
receivers as observed from the secondary transmitters. It is worth emphasizing that the structure of the 
null constraints in (|3j) is much more general than the three cases mentioned above, as it can incorporate 
any combination of the frequency, time and space coordinates. 

The use of the spatial domain can greatly improve the capabilities of cognitive users, as it allows 
them to transmit over the same frequency band but without interfering. This is possible if the secondary 
transmitters have an antenna array and use a beamforming that puts nulls over the directions identifying 
the primary receivers. Of course, this requires the identification of the primary receivers, a task that is 
much more demanding than the detection of primary transmitters [4J. As an example, there are some 
recent works showing that, in the application of CR over the spectrum allocated to commercial TV, 
one might exploit the local oscillator leakage power emitted by the RF front end of the TV receiver to 
locate the receivers [18]. Of course, in such a case, the detection range is quite short and this calls for 
a deployment of sensors very close to the potential receivers. A different scenario pertains to cellular 
systems. In such a case, the mobile users might be rather hard to locate and track. However, the base 
stations are relatively easier to identify. Hence, in a cellular system operating in a time-division duplexing 
(TDD) mode, the secondary users could exploit the time slot allocated for the uplink channel and put a 
null in the direction of the base stations. This would avoid any interference towards the cellular system 
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users, without the need of tracking the mobile users. 

The soft shaping constraints expressed in (|4]) and j5j) represent a constraint on the total average and 
peak average power radiated (projected) along the directions spanned by the column space of matrix 
G q . They are a relaxed form of jHJ) and can be used to keep the portion of the interference temperature 
generated by each secondary user q under the desired value. In fact, under flU)-©, the secondary users 
are allowed to transmit over some subchannels occupied by the primary users, but only provided that 
the interference that they generate falls below a prescribed threshold. For example, in a MIMO setup, 
the matrix G q in (J4j) would contain, in its columns, the steering vectors identifying the directions of the 
primary receivers. 

Within the assumptions made above, invoking the capacity expression for the single user Gaussian 
MIMO channel— achievable using random Gaussian codes by all the users— the maximum information rate 
on link q for a given set of users' covariance matrices Qi, . . . , Qq, is [19] 

R q {Q q , Q-q) = log det (I + Hg q R,ZlHqqQq) (6) 

where 

R-g — Rn.q + y^H rg Q r H,^ (7) 

r^q 

is the MUI plus noise covariance matrix observed by user q and Q q = (Q r ) r ^ q is the set of all the users' 
covariance matrices, except the q-th one. Observe that R_ g depends on the strategies Q- q of the other 
players. 

3 Resource Sharing among Secondary Users based on Game Theory 

Given the multiuser nature of the scenario described above, the design of the optimal transmission strate- 
gies of secondary users would require a multiobjective formulation of the optimization problem, as the 
information rate achieved on each secondary user's link constitutes a different single objective func- 
tion. The globally optimal solutions of such a problem— the Pareto optimal surface of the multiobjective 
problem— would define the largest rate region achievable by secondary users, given the power constraints 
Co.l-Co.4: the rate vector profile R(Q*) = [i?i(Q*) 5 . . . ,Rq(Q*)] is Pareto optimal if there exists no 
other rate profile R(Q) that dominates R(Q*) component-wise, i.e., R(Q*) > R(Q), for all feasible Q's, 
where at least one inequality is strict. 

Unfortunately, the computation of the rate region is analytically intractable and thus not applicable 
in a cognitive radio scenario, since every scalar /multiobjective optimization problem involving the rates of 
secondary users in j6]) is not convex (implied from the fact that the rates R q (Q) are nonconcave functions 
of the covariance matrices Q). Furthermore, even in the simpler case of transmissions over SISO parallel 
channels, the network utility maximization (NUM) problem based on the rates functions J6)) has been 
proved in [24j to be a strongly NP-hard problem, under various practical settings as well as different 
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choices of the system utility function (e.g., sum-rate, weighted sum-rate, geometric rate-mean). Roughly 
speaking, this means that there is no hope to obtain an algorithm, even centralized, that can efficiently 
compute the exact globally optimal solution. Although in theory, the rate region could be still found by an 
exhaustive search through all possible feasible covariance matrices, the computational complexity of this 
approach is prohibitively high, given the large number of variables and users involved in the optimization. 
The situation is particularly critical in CR systems, where the cognitive users sense a very large spectrum. 
Consequently, suboptimal algorithms have been proposed in the literature to solve special cases of the 
proposed optimization [20]- [23], most of them dealing with the maximization of the (weighted) sum-rate 
in SISO frequency-selective interference channels (obtained from our general model when the channel 
matrices are diagonal, the covariance matrices reduce to the power allocation vectors, and the null/soft 
shaping constraints are removed) [20j[21]. Due to the nonconvex nature of the problem, these algorithms 
either lack global convergence or may converge to poor spectrum sharing strategies. 

Furthermore, even if one decides to employ a suboptimal method, e.g., [20]-[23j. the algorithms are 
not suitable for CR systems as they are centralized and thus cannot be implemented in a distributed way. 
These techniques require a central authority (or node in the network) with knowledge of the (direct and 
cross-) channels to compute all the transmission strategies for the different nodes and then to broadcast 
the solution. This scheme would clearly pose a serious implementation problem in terms of scalability of 
the network and amount of signaling to be exchanged among the nodes, which makes such an approach 
not appealing in the scenario considered in this paper. 

To overcome the above difficulties and reach a better trade-off between performance and complexity, we 
shift our focus to a different notion of optimality: the competitive optimality criterion; which motivates a 
game theoretical formulation of the system design. Using the concept of NE as the competitive optimality 
criterion, the resource allocation problem among secondary users is then cast as a strategic noncooperative 
game, in which the players are the secondary users and the payoff functions are the information rates on 
each link: Each secondary user q competes against the others by choosing the transmit covariance matrix 
Q 9 (i.e., his strategy) that maximizes his own information rate .R 9 (Q 9 ,Q_ 9 ) in ([6]), given constraints 
imposed by the presence of the primary users, besides the usual constraint on transmit power. A NE of 
the game is reached when each user, given the strategy profiles of the others, does not get any rate increase 
by unilaterally changing his own strategy. The first question to answer under such framework is whether 
such an overall dynamical system can eventually converge to an equilibrium point, while preserving the 
QoS of primary users. The second basic issue is if the optimal strategies to be adopted by each user can 
be computed in a totally decentralized way. We address both questions in the forthcoming sections. 

For the sake of simplicity, we start considering only constraints Co.l and Co. 2. These constraints 
are suitable to model interweave communications among secondary users where, in general, there are 
restrictions on when and where they may transmit (this can be done using the null constraints Co. 2). 
Then, we allow underlay and interweave communications simultaneously, by including in the optimization 
also interference constraints Co. 3 and Co. 4. 
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3.1 Rate maximization game with null constraints 

Given the rate functions in j6|) and constraints Co.l-Co.2, the rate maximization game is formally defined 



as: 



maximize i?„(Q a ,Q_„) 
(Ql): Vq=l,---,Q, (8) 

subject to Tr (Q ? ) < P q , Uf Q q = 

where Q is the number of players (the secondary users) and R q (Q q , Q-q) is the payoff function of player 

q, defined in j6]). Without the null constraints, the solution of each optimization problem in J8]) would 

lead to the well-known MIMO waternlling solution [19]. The presence of the null constraints modifies the 

problem and the solution for each user is not necessarily a waternlling anymore. Nevertheless, we show 

now that introducing a proper projection matrix the solutions of J8]) can still be efficiently computed via 

a waterfilling-like expression. To this end, we rewrite game Q\ in a more convenient form as detailed next. 

Introducing the projection matrix P-^(u q )- L = I — U^U^Uq) -1 !^ (the orthogonal projection onto 

TZ(\J q )' L , where TZ(-) is the range space operator), it follows from the constraint XJ^Q q = that any 

optimal Q g in JBJ) will always satisfy: 

Qg = P 7e(u,)- L Qg P 7^(u,)- L - ( 9 ) 
The game Q\ can then be equivalently rewritten as: 

maximize log det ( I + H^RI g Hq g C 
subject to Tr(Q 9 )<P g Vg = l,- ••,<?, (10) 

Qg = p ^(u g )-LQg p -R.(u 9 )^ 

where each H rg = H rg P^ Ur )± is a modified channel and R-g = R„ 9 + ^ H rg Q r H,^. At this point, 

r^q 

the problem can be further simplified by noting that the constraint Q g = P^(u^) Qg P 7£(u^) i n (fTOjl is 
redundant. The final formulation then becomes: 

maximize log det ( I + H^RIgH„gQ„ ) 

V 1 Vg = l,--- ,Q. (11) 

subject to Tr(Qg) < P q 

This is due to the fact that, for any user q, any optimal solution Q* in (fTTjl — the MIMO waternlling solution 
[13]— will be orthogonal to the null space of H gg , whatever R_ g is, implying Q* = P7?.(u 9 )- L QgP7e(u g ) ± • 
Building on the equivalence of (jSJ) and (fTTI) . we can apply the results in [13J to the game in (fTTj) and 
derive the structure of the Nash equilibria of game Gi, as detailed next. 

Nash equilibria of game Q\\ Game Q\ always admits a NE, for any set of channel matrices, transmit 
power of the users, and null constraints, since it is a concave game (the payoff of each player is a concave 
function in his own strategy and each admissible strategy set is convex and compact) [13]. Moreover, 
it follows from (PTT1) that all the Nash equilibria of Q\ satisfy the following set of nonlinear matrix- value 
fixed-point equations [13j: 

Ql = WF, (Hf 9 R: g (Q* ? )Hgg) 4 Diag (p£) Wf , Vg = 1, • • • , Q, (12) 

9 



where we made explicit the dependence of R_ g on Q- q as R_ 9 (Q* ); the W* = W 9 (Q* 9 ) is the 
semi-unitary matrix with columns equal to the eigenvectors of matrix H^RI*(Q* q )fl qq corresponding 
to the positive eigenvalues X* k = A (?j fe(Q* (J ), with R_ g (Q_ g ) defined in J7)); and the power allocation 
Pq = Pg(Q-g) satisfies the following simultaneous waterfilling equation: for all k and q, 

P* q (k) = (^ q ~ , (13) 

with (x) + = max(0, x) and \x q chosen to satisfy the power constraint Y^kPqi^) = 

Interestingly, the solution (fl~2j) shows that the null constraints in the transmissions of secondary users 
can be handled without affecting the computational complexity: The optimal transmission strategy of 
each user q can be efficiently computed via a MIMO waterfilling solution, provided that the original 
channel matrix H gq is replaced by H gg . 

This result has an intuitive interpretation: To guarantee that each user q does not transmit over a 
given subspace (spanned by the columns of U g ), whichever the strategies of the other users are, while 
maximizing his information rate, one only needs to induce in the channel matrix H qq a null space that 
coincides with the subspace where the transmission is not allowed. This is precisely what is done by 
introducing the modified channel H ? q- 

The waterfilling-like structure of the Nash equilibria as given in (fT2l) along with the interpretation 
of the MIMO watefilling solution as a matrix projection onto a proper convex set as given in [13| play 
a key role in studying the uniqueness of the NE and in deriving conditions for the convergence of the 
distributed algorithms described in Section [H The analysis of the uniqueness of the NE goes beyond the 
scope of this paper and it is addressed in [14]. What is important to remark here is that, as expected, 
the conditions guaranteeing the uniqueness of the NE impose a constraint on the maximum level of MUI 
generated by secondary users that may be tolerated in the network. But, interestingly, the uniqueness of 
the equilibrium is not affected by the interference generated by the primary users. 

3.2 Rate maximization game with null constraints via virtual noise shaping 

In this section, we show that an alternative approach to impose null constraints Co. 2 on the transmissions 
of secondary users passes through the introduction of virtual interferers. The idea behind this alterna- 
tive approach can be easily understood if one considers the transmission over SISO frequency-selective 
channels, where all the channel matrices have the same eigenvectors (the FFT vectors): to avoid the use 
of a given subchannel, it is sufficient to introduce a "virtual" noise with sufficiently high power over that 
subchannel. The same idea cannot be directly applied to the MIMO case, as arbitrary MIMO channel 
matrices have different right/left singular vectors from each other. Nevertheless, we show how to design 
the covariance matrix of the virtual noise (to be added to the noise covariance matrix of each secondary 
receiver), so that the all the Nash equilibria of the game satisfy the null constraint Co. 2 along the specified 
directions. 
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Let us consider the following strategic noncooperative game: 

maximize log det (i + H^.R~ „H 00 Q ) 
{9a) i V 9 = l > -,0, (14) 

subject to Tr (Q g ) < P q 

where 

R_,, a 4 R + aU.Uf = R„, + H r? Q r H^ + aU 9 Uf , (15) 

r^g 

denotes the MUI-plus-noise covariance matrix observed by secondary user q, plus the covariance matrix 
aXJqXJ^ of the virtual interference along TZ(JJ q ), where U g is a tall matrix and a is a positive constant. 
Our interest is on deriving the asymptotic properties of the solutions of Q a , as a — > +oo. To this end, 
we introduce the following intermediate definitions first. For each q, define the tall matrix U^- such that 
7l(XJq) = IZiXJq)- 1 , and the modified channel matrices 

H rq = tifH rq Vr, q=l,--- ,Q- (16) 

We then introduce the auxiliary game Goo, defined as: 

maximize log det ( I + H^RI*H 00 Q ) 
(Goo): V ' Vg = l,---,Q > (17) 

subject to Tr(Qq) < P q 

where 

R-g — U^R n<7 U^ + H rg Q r H^. (18) 

It can be shown that games G a and Goo are asymptotically equivalent in the sense specified next. 

Nash equilibria of games G a and Goo- Games Ga and Goo always admit a NE, for any set of channel 
matrices, power constraints, and a > 0. Moreover, under mild conditions guaranteeing the uniqueness of 
the NE of both games (denoted by Q* and QJo, respectively), we have: 

lim Q* = Q^, (19) 

a— >oo 

i.e., the NE of Ga asymptotically coincides with that of Goo- 

Observe that, similarly to game Gi, also in games Ga and Goo, the best-response of each player can be 
efficiently computed via MIMO waterfilling-like solutions, and the Nash equilibria of both games satisfy 
a simultaneous waterfilling equation. 

Using (fl~9j) . one can derive the asymptotic properties of the (unique) NE of game Ga as a — > oo, 
through the properties of the equilibrium of Goo- Following a similar approach as in Section [3TTT one 
can show that each Qg jOC satisfies the following condition 

UfQ* iOO = 0, with U^H- 1 ^. (20) 

Condition (|20j) provides, for each user q, the desired relationship between the directions of the virtual 
noise to be introduced in the noise covariance matrix of the user (see ([18]) ) — the matrix Ug — and the real 



11 



directions along with user q will not allocate any power, i.e., the matrix U ? . It turns out that if user q is 
not allowed to allocate power along U q , it is sufficient to choose in (fT8|) JJ q = H qq U q . 

Since the existence and uniqueness of the NE of game Q a do not depend on a, the (unique) NE of 
Q a (that in general will depend on the value of a) can be reached using the asynchronous algorithms 
described in Section HI irrespective of the value of a. Thus, for sufficiently large values of a, the NE of 
Q a tends to satisfy condition (120]) : which provides an alternative way to impose constraint Co. 2. 

3.3 Rate maximization game with soft and null constraints 

We focus now on the rate maximization in the presence of both null and soft shaping constraints. The 
resulting game can be formulated as follows: 

maximize R q (Q q , Q- q ) 
. subject to Tr (Gf Q,G 9 ) < Pf e V g = 1, • • • , Q. (21) 

Amax (G?QgG„) < pP° av 



U^Q 9 = 



r <? ^q^qj -i 1 q 

'q ^q 

We assume w.l.o.g. that each G q is a full-row rank matrix, so that the soft shaping constraint in ((21]) 
imposes a constraint on the average transmit power radiated by user q in the whole space. 

The soft constraints in (|2Tj) are the result of a constraint on the overall interference temperature limit 
imposed by the primary users [2]. Typically, the most stringent conditions between the power constraints 
Co.l and Co. 3 is the soft shaping constraint Co. 3. This motivates the absence in ((21]) of the power 
constraint Co.l, although it could also be considered. 

Nash equilibria of game 0%: We can derive the structure of the Nash equilibria of game Q2, similarly 
to what we did for game Q\. For each q G f2, define the tall matrix U,j = G q XJ q , where G| denotes the 
Monroe-Penrose pseudoinverse of G q [25], introduce the projection matrix P^u,,^ = I — U g (U^U g ) _1 
(the orthogonal projection onto ^(Ug)^) and the modified channel matrices 

H rg = H rg Gj. P-^Ur^ ' r > 9 = 1) ' ' ' ; Q- (22) 

Using the above definition, we can now characterize the Nash equilibria of game G2, as shown next. 

The game Q2 admits a NE, for any set of channel matrices and null/soft shaping constraints. Moreover, 
every NE satisfies the following set of nonlinear matrix-value fixed-point equations: 

Q- = GrWF,(H>:; ( Q % )H M )G! 



Gl V* diag (p*) V* H G 



where V* = V g (Ql 9 ) is the semi-unitary matrix with columns equal to the eigenvectors of matrix 
H^RIg(Q* q )H qq , with R_ g (Q_ g ) defined in (jTj) , corresponding to the L q = rank(H gg ) positive eigen- 
values X q k = Aq^Q* q ), and the power allocation p* = p 9 (Q* g ) satisfies the following simultaneous 
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waterfilling equation: for all k and q, 

^pcak 



p*Jk) 



if P$ eak L n > P? ve 



pP eak ) otherwise, 



o " (24) 



ppcak 

where [-] 9 denotes the Euclidean projection onto the interval [0, Pq ea ] and \i q is chosen to satisfy the 
power constraint Y2kPq(^) = ^"g™ ( see i e -S-i [23] f° r practical algorithms to compute such a /J, q ). 

The structure of the Nash equilibria in ([23]) states that the optimal transmission strategy of each user 
leads to a diagonalizing transmission with a proper power allocation, after pre/post multiplication of the 
waterfilling solution by matrix G|. Similarly to Gi, the conditions for the uniqueness of the NE of game 
Q2 can be obtained, building on the interpretation of the waterfilling solutions in (|23l) as matrix projection 
[T3] . As expected, the NE of the game is unique, provided that the interference generated by secondary 
users is not too high. 



4 MIMO Asynchronous Iterative Waterfilling Algorithm 

In Section [3] we have shown that the optimal resource allocation among secondary users in hierarchical 
cognitive networks corresponds to an equilibrium of the system, where all the users have maximized 
their own rates, without hampering the communications of primary users. Since there is no reason to 
expect a system to be initially at the equilibrium, the fundamental problem becomes to find a procedure 
that reaches such an equilibrium from non-equilibrium states. In this section, we focus on algorithms 
that converge to these equilibria. Since we are interested in a decentralized implementation, where no 
signaling among secondary and primary users is allowed, we consider only totally distributed iterative 
algorithms, where each user acts independently of the others to optimize his own transmission strategy 
while perceiving the other active users as interference 

More specifically, to reach the Nash equilibria of the games introduced in the previous section, we 
propose a fairly general distributed and asynchronous iterative algorithm, called asynchronous Iterative 
WaterFilling Algorithm (IWFA). In this algorithm, all secondary users maximize their own rate [via the 
single user MIMO waterfilling solution (fT2j) for game Qi, ([23]) for game Q2, and the classical MIMO 
waterfilling solution for games Q a and Qoo\ in a totally asynchronous way, while keeping the temperature 
noise levels in the licensed bands under the required threshold [2]. According to the asynchronous updating 
schedule, some users are allowed to update their strategy more frequently than the others, and they might 
even perform these updates using outdated information on the interference caused by the others. 

Before introducing the proposed asynchronous MIMO IWFA, we need the following preliminary defi- 
nitions. We assume, without loss of generality, that the set of times at which one or more users update 
their strategies is the discrete set T = N+ = {0, 1,2,...}. Let Qq denote the covariance matrix of the 
vector signal transmitted by user q at the n-th iteration, and let T q C T denote the set of times n at 
which is updated (thus, at time n T q , is left unchanged). Let r^(n) denote the most recent 
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time at which the interference from user r is perceived by user q at the n-th iteration (observe that rf(n) 
satisfies < T$(n) < n). Hence, if user q updates his own covariance matrix at the n-th iteration, then he 
chooses his optimal Q q n \ according to (fl~2|) for game Q\ and ([23]) for game G2, and using the interference 
level caused by 

Q(-r q {n)) A 

Some standard conditions in asynchronous convergence theory that are fulfilled in any practical imple- 
mentation need to be satisfied by the schedule {T~r(n)} and {Tq}; we refer to [13] for the details. Using 
the above notation, the asynchronous MIMO I WFA is formally described in Algorithm 1 below, where 
the mapping in (|27|) is defined as 

T,(Q_ 9 ) = WF g (h*RiJH w ) , q=l,-..,Q, (26) 

with WFq (•) given in (fT2j) if the algorithm is applied to game Q\, and it is defined as 

T 9 (Q_ 9 ) 4 G f WF, (Hf.RljH^) G», q = 1, ■ ■ ■ , Q, 

with WF 9 (•) given in ([23]) if the algorithm is applied to game Q?,- The mapping T 9 (Q_ 9 ) reduces to the 
classical MIMO waterfilling solution [19] if games Q a and Goo are considered. 



Qi 



(r?(n)) 



9-1 



Q 



(25) 



Algorithm 1: MIMO Asynchronous I WFA 

Set n = and Qq ^ = any feasible point; 
f or n = : N it 

1 Q? n \ otherwise; 

end 



Convergence of the asynchronous I WFA is studied in [T3], [H] (see also [IU [TJ] for special cases of the 
algorithm), where it was proved that the algorithm converges to the NE of the proposed games under 
the same conditions guaranteeing the uniqueness of the equilibrium. The proposed asynchronous I WFA 
contains as special cases a plethora of algorithms, each one obtained by a possible choice of the schedule 
{r?(n)}, {T q }. The sequential [21 QU EH EH] and simult aneous [H]-[l3] IWFAs are just two examples 
of the proposed general framework. The important result stated in [TT]-[l3] is that all the algorithms 
resulting as special cases of the asynchronous MIMO IWFA are guaranteed to reach the unique NE of 
game under the same set of convergence conditions, since convergence conditions do not depend on the 
particular choice of {T q } and {rr(n)} [13] . 

Moreover all the algorithms obtained from Algorithm 1 have the following desired properties: 

- Low complexity and distributed nature: Even in the presence of null and/or shaping constraints, 
the best response of each user q can be efficiently and locally computed using a MIMO waterfilling based 
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Figure 3: Simultaneous vs. sequential IWFA: rates of secondary users versus iterations, obtained by the sequential 
IWFA (dashed- line curves) and simultaneous IWFA (solid- line curves). 

solution, provided that each channel H qq is replaced by the modified channel H gg (if game Q\ is considered) 
or H qq (if game Q2 is considered). Thus, Algorithm 1 can be implemented in a distributed way, since each 
user only needs to measure the overall interference-plus-noise covariance matrix R_<j and waterfill over 
H^RIgHqg [or over H^RljH^]. 

- Robustness: Algorithm 1 is robust against missing or outdated updates of secondary users. This 
feature strongly relaxes the constraints on the synchronization of the users' updates with respect to those 
imposed, for example, by the simultaneous or sequential updating schemes [TlJ-[13j. 

- Fast convergence behavior: The simultaneous version of the proposed algorithm converges in a very 
few iterations, even in networks with many active secondary users. As an example, in Figure [3] we show 
the rate evolution of the of 3 links out 8 secondary users, corresponding to the sequential IWFA and 
simultaneous IWFA as a function of the iteration index. As expected, the sequential IWFA is slower 
than the simultaneous IWFA, especially if the number of active secondary users is large, since each user 
is forced to wait for all the users scheduled in advance, before updating his own covariance matrix. This 
intuition is formalized in [llj, where the authors provided the expression of the asymptotic convergent 
factor of both the sequential and simultaneous IWFAs. 

- Control of the radiated interference: Thanks to the game theoretical formulation including null 
and/or soft shaping constraints, the proposed asynchronous IWFA does not suffer of the main drawback 
of the classical sequential IWFA [27], i.e., the violation of the interference temperature limits [2]. 

Figure 0] shows an example of the optimal resource allocation based on the game theoretical for- 
mulation Qi, for a cognitive MIMO network composed by two primary users and two secondary users, 
sharing the same spectrum and space. Secondary users are equipped with four transmit/receive an- 
tennas, placed in uniform linear arrays critically spaced at half of the wavelength of the passband 
transmitted signal. For the sake of simplicity, we assumed that the channels between the transmitter 
and the receiver of the secondary users have three physical paths (one line-of-sight and two reflected 
paths) as shown in Figure HJa). To preserve the QoS of primary users' transmissions, null constraints 
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Figure 4: Optimal transmit beamforming patterns at the NE of game Gi [subplot (b)] for a cognitive MIMO 
network composed by two primary and two secondary users [subplot (a)]. 

are imposed to secondary users in the (line-of-sight) "directions" of primary users' receivers [see sub- 
plot (a)]. For the scenario shown in the figure, one null constraint for each player is imposed along 
the transmit directions = ir/2 and </> 2 = — 57r/12. This can be done choosing for each player q 
the matrix XJ q in (j2"T1) coinciding with the spatial signature vector in the transmit direction (f) q , i.e., 
Ug = [l,exp(-j27rA t(j sin(0 9 )),exp(-j27r2A t(3 sin(^)),exp(-j27r3A t9 sin((/) g ))] T , with A tq = 1/2 denot- 
ing the normalized (by the signal wavelength) transmit antenna separation and q = 1,2. In Figure E](b) , 
we plot the transmit beamforming patterns, associated to the eigenvectors of the optimal covariance ma- 
trix of the two secondary users at the NE, obtained using Algorithm 1. In each radiation diagram plot, 
solid (blue) and dashed (black) line curves refer to the two eigenvectors corresponding to the nonzero 
eigenvalues (arranged in increasing order) of the optimal covariance matrix [recall that, because of the 
null constraints, the equivalent channel matrix H qq in (|2"T1) has rank equal to 2]. Observe that the null 
constraints guarantee that at the NE no power is radiated by the two secondary transmitters along the 
directions 4> 1 (for transmitted one) and cf> 2 (for transmitted two), showing that in the MIMO case, the or- 
thogonality among primary and secondary users can be reached in the space rather than in the frequency 
domain, implying that primary and secondary users may share frequency bands, if this is allowed by FCC 
spectrum policies. 

5 Special Cases 

The MIMO game theoretic formulation proposed in the previous sections provides a general and unified 
framework for studying the resource allocation problem based on rate maximization in hierarchical CR 
networks, where primary and secondary users coexist. In this section, we specialize the results to two 
scenarios of interest: 1) the spectrum sharing problem among primary and secondary users transmitting 
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over SISO frequency-selective channels; and 2) the MIMO transceivers design of heterogeneous systems 
sharing the same spectrum over unlicensed bands. 

5.1 Spectrum sharing over SISO frequency-selective channels with spectral mask 
constraints 

The block transmission over SISO frequency-selective channels is obtained from the I/O model in ([[]), 
when each channel matrix H rq is a N x N Toeplitz circulant matrix, R„ q is a N x N diagonal matrix N is 
the length of the transmitted block (see, e.g., [ID])- This leads to the following eigendecomposition for each 
channel H rq = WB rq W H , where W is the normalized IFFT matrix, i.e., [W] y = ef^-W-W 1 * /y/N 
for i, j = 1, . . . , N and D rq is & N x N diagonal matrix, where [D rg ] fcfc = H rq (k) is the frequency-response 
of the channel between source r and destination q. Within this setup, we focus on game Q\ given in (jSJ), 
but similar results could be obtained if game Q2, Qa or Qoo were considered instead. In the case of SISO 
frequency-selective channels, game Q\ can be rewritten as: 

maximize logdet (I + H^R~ H ?9 Q g ) 

subject to Tr(Q,)<P, Vg = l,---,Q, (28) 

[W H Q,W] tt <p^(fc), Vk = l r --,N, 

where {p™ ax (&0} is the set of spectral mask constraints, that can be used to impose shaping (and thus also 
null) constraints on the transmit power spectral density (PSD) of secondary users over licensed/unlicensed 
bands. 

Nash equilibria: The solutions of the game in (|28j) have the following structure [10J: 

Q£ = WDiag(p*)W H , Wq=l,---,Q, (29) 
where p* = (p* (k))^ =1 is the solution to the following set of fixed-point equations 

p* = wf 9 (p* q ) , Vg=l,-.. ,Q, (30) 

with the waterfilling vector operator wf g (•) defined as 

, k = l,---,N, (31) 



where (i q is chosen to satisfy the power constraint with equality ^2 k p q {k) = P q . 

Equation (|29l) states that, in the case of SISO frequency-selective channels, a NE is reached using, 
for each user, a multicarrier strategy (i.e., the diagonal transmission strategy through the frequency 
bins), with a proper power allocation. This simplification with respect to the general MIMO 
consequence of the property that all channel Toeplitz circulant matrices are diagonalized by the same 
matrix, i.e., the IFFT matrix W, that does not depend on the channel realization. 



[wf g (p_ 



i + E 



r9 £q\"rq(k)\ Pr(k) 



\H qq {k)Y 
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V max 1 _ J \ < 1, Wr = l,---,Q, and VA: = 1, • • • ,N, (33) 

^ * Lff M (*:) 2 te. - - 



Interestingly, multicarrier transmission with a proper power allocation for each user is still the opti- 
mal transmission strategy if in (I28p instead of the information rate, one considers the maximization of 
the transmission rate using finite order constellations and under the same constraints as in (l28l) plus a 
constraint on the average error probability. Using the gap approximation analysis, the optimal power 
allocation is still given by the waterfilling solution ([31]) . where each channel transfer function \H qq (k)\ 2 is 
replaced by \H qq (k)\ 2 /T q , where T g > 1 is the gap [10] , The gap depends only on the family of constella- 
tion and on error probability constraint P e ,q', for M-QAM constellations, for example, the resulting gap 
is = (Q- 1 (P e , g /4)) 2 /3 (see, e.g., [29]). 

Reaching a NE of the game in (|28l) satisfies a competitive optimality principle, but, in general, multiple 
equilibria may exist, so that one is never sure about which equilibrium is really reached. Sufficient 
conditions on the MUI that guarantee the uniqueness of the equilibrium have been proposed in the 
literature [I0]-[I2] and [27J [28]. Among all, one of the two following conditions is sufficient for the 
uniqueness of the NE: 

\H (k)\ 2 d 2 

Vmax 1 rqK % < 1, Mq= 1, • • • , Q, and V& = 1, • • • , N, (32) 
^ * \H qq (k)\ 2 d 2 rq 

|i?j.q(A;)| d qq 

r+q " \H qq (k)\ 2 ^ 

where we have introduced the normalized channel transfer functions H rq (k) = H rq (k)/d 2 , Mr, q, with d rq 
indicating the distance between transmitter of the r-th link and the receiver of the q-th link. From (1321) - 
(1331) . it follows that, as expected, the uniqueness of NE is ensured if secondary users are sufficiently far 
apart from each other. In fact, from (l3"2"l) - (l3"3[) for example, one infers that there exists a minimum distance 
beyond which the uniqueness of NE is guaranteed, corresponding to the maximum level of interference 
that may be tolerated by the users. Specifically, condition (13"2[) imposes a constraint on the maximum 
amount of interference that each receiver can tolerate; whereas (I33p introduces an upper bound on the 
maximum level of interference that each transmitter is allowed to generate. Interestingly, the uniqueness 
of the equilibrium does not depend on the interference generated by the transmissions of primary users. 

Asynchronous IWFA: To reach the equilibrium of the game, secondary users can perform the asyn- 
chronous IWFA based on the mapping in ([HJ. This algorithm can be obtained directly from Algorithm 
1, as special case. It was proved in [12] that, e.g., under conditions ([32]) - ([33]) . the asynchronous IWFA 
based on mapping (|3Tj) converges to the unique NE of game in (l28l) as Nit— ► oo, for any set of feasible 
initial conditions and updating schedule. 

In Figure [5j we show an example of the optimal power allocation in SISO frequency-selective channels 
at the NE, obtained using the proposed asynchronous IWFA, for a CR system composed by one primary 
user [subplot (a)] and two secondary users [subplot (b)], subject to null constraints over licensed bands, 
spectral mask constraints and transmit power constraints. In each plot, solid and dashed-dot line curves 
refer to optimal PSD of each link and PSD of the MUI plus thermal noise, normalized by the channel 
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Figure 5: Spectrum sharing among one primary [subplot (a)] and two secondary users [subplot (b)]: Optimal 
PSD of each link (solid lines), and PSD of the MUI-plus-thermal noise normalized by the channel transfer function 
square modulus of the link (dashed-dot line) . 

transfer function square modulus of the link, respectively. In this example, there is a band A (from 50 
to 300 frequency bins) allocated to an active primary user; there is then a band B (from 300 to 400 
frequency bins) allocated to licensed users, but temporarily unused; the rest of the spectrum, denoted as 
C, is vacant. The temporarily void band B can be utilized by secondary users, provided that they do not 
overcome a maximum tolerable spectral density. The optimal power allocations shown in Figure [5] are 
the result of running the simultaneous I WFA. We can observe that the secondary users do not transmit 
over band A and they allocate their power over both bands B and C, respecting a power spectral density 
limitation over band B. 



5.2 MIMO transceiver design of heterogeneous systems in unlicensed bands 

We consider now on a scenario where multiple unlicensed MIMO cognitive users share the same unlicensed 
spectrum and geographical area. The availability of MIMO transceivers clearly enriches the possibilities 
for spectrum sharing as it adds the extra spatial degrees of freedom. In unlicensed bands, there are no 
interference constraints to be satisfied by the users. Thus, the game theoretical formulation as given in 
(|Sj|, without considering the null constraints, seems the most appropriate to study the resource allocation 
problem in this scenario. In the following we refer to game Q\ assuming tacitly that the null constraints 
are removed. 

Similarly to the SISO case, sufficient conditions for the uniqueness of the NE are given by one of the 
two following set of conditions (more general conditions are given in [13J): 

Low MUI received: ^ p (H^H^E^/H^) < 1, Vg = 1, • • • ,Q, (34) 

r^q 

Low MUI generated: ^^(H^H^H^H^) < 1, Vr = 1, ■ • • ,Q. (35) 

q+r 

Conditions (|3"4"1) - ([35"1) quantify how much MUI can be tolerated by the systems to guarantee the uniqueness 
of the NE. Interestingly, (j32j) - ([33"1) and most of the conditions known in the literature [TTJ[27l[28] for the 
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uniqueness of the NE of the rate-maximization game in SISO frequency-selective interference channels 
and OFDM transmission come naturally from (|34l) - ([35l) as special cases. 

The Nash equilibria of game Q\ can be reached using the asynchronous IWFA described in Algorithm 1, 
whose convergence is guaranteed under conditions (|34l) - ([35j) . for any set of initial conditions and updating 
schedule of the users. 




Figure 6: Sum-Rate of the users versus the inter-pair distance d rq /d qq ; d rq = d qr , d rr = d qq = 1, Vr, q, for different 
numbers of transmit /receive antennas. 

In Figure [6] we show an example of the benefits of MIMO transceivers in the cognitive radio context. 
We plot in the figure the sum-rate of a two-user frequency-selective MIMO system as a function of the 
inter-pair distance among the links, for different number of transmit/receive antennas. The rate curves are 
averaged over 500 independent channel realizations, whose taps are simulated as i.i.d. Gaussian random 
variables with zero mean and unit variance. For the sake of simplicity, the system is assumed to be 
symmetric, i.e., the transmitters have the same power budget and the interference links are at the same 
distance (i.e., d rq = d qr , Vg,r), so that the cross channel gains are comparable in average sense. From the 
figure one infers that, as for isolated single-user systems or multiple access/broadcast channels, also in 
MIMO interference channels, increasing the number of antennas at both the transmitter and the receiver 
side leads to a better performance. The interesting result, coming from Figure EJ is that the incremental 
gain due to the use of multiple transmit/receive antennas is almost independent of the interference level 
in the system, since the MIMO (incremental) gains in the high-interference case (small values of d rq /d qq ) 
almost coincide with the corresponding (incremental) gains obtained in the low-interference case (large 
values of d rq /d qq ), at least for the system simulated in Figure [H This desired property is due to the 
fact that the MIMO channel provides more degrees of freedom for each user than those available in the 
SISO channel, that can be explored to find out the best partition of the available resources for each user, 
possibly cancelling the MUI. 
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6 Conclusion and Directions for Further Developments 



In this paper we have proposed a signal processing approach to the design of CR systems, using a 
competitive optimality principle, based on game theory. We have addressed and solved some of the 
challenging issues in CR, namely: 1) the establishment of conditions guaranteeing that the dynamical 
interaction among cognitive nodes, under constraints on the transmit spectral mask and on interference 
induced to primary users, admits a (possibly unique) equilibrium; and 2) the design of decentralized 
algorithms able to reach the equilibrium points, with minimal coordination among the nodes. We have 
seen how basic signal processing tools such as subspace projectors play a fundamental role. The spectral 
mask constraints have been in fact used in a very broad sense, meaning that the projection of the 
transmitted signal along prescribed subspaces should be null (null constraints) or below a given threshold 
(soft constraints). The conventional spectral mask constraints can be seen as a simple case of this general 
set-up, valid for SISO channels and using as subspaces the space spanned by the IFFT vectors with 
frequencies falling in the guard bands. This general setup encompasses multiantenna MIMO systems, 
which is particularly useful for CR, as it provides the additional spatial degrees of freedom to control the 
interference generated by the cognitive users. 

Of course, this field of research is full of interesting further directions worth of investigation. The 
NE points derived in this paper were dictated by the need of finding totally decentralized algorithms 
with minimal coordination among the nodes. However, the NE points may not be Pareto-efficient. This 
raises the issue of how to move from the NE towards the Pareto optimal trade-off surface, still using a 
decentralized approach. Game theory itself provides a series of strategies to move from inefficient Nash 
equilibria towards Pareto-efficient solutions, still using a decentralized approach, through, for example, 
repeated games, where the players learn from their past choices [9]. Examples of such games are the auction 
games, where the auctioneer (primary users) dynamically determine resource allocation and prices for the 
bidders (secondary users), depending on traffic demands, QoS and supply/demand curves, as evidenced 
in a series of works (see, e.g., [30j [HJ E2]). Repeated games may also take the form of negotiations 
between primary and secondary users, with primary users willing to lease part of their spectrum to 
secondary users, under suitable remunerations [16] or under the availability given by secondary users 
to establish cooperative links with the primary users to improve their QoS [33]. Competitive pricing for 
spectrum sharing was also proposed as an oligopoly market where a few primary users offer spectrum access 
opportunities to secondary users [34]. An interesting issue will be the integration of our asynchronous 
IWFA in repeated (auction) games, where the optimization considers a set of primary users offering the 
lease of portion of their resources to a set of secondary users, as a function of traffic demands, QoS 
requirements and physical constraints. 

Our search for the uniqueness conditions of the NE and the convergence conditions of our proposed 
algorithms forced us to simplify the model. For example, we assumed that each receiver has an error-free 
short-term prediction of the channel. This assumption was necessary for the mathematical tractability 
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of the problem and to be able to provide closed-form expressions of our findings. This is useful to gain a 
full understanding of the problem, without relying on simulation results only. However, in practice, the 
transmitter is only able to acquire an estimate affected by errors and, based on that, to form a prediction 
of the short term future evolution. An interesting extension of the presented approach consists then in 
taking into account the effects of estimation errors and developing robust strategies. This is particularly 
relevant in CR systems because the strategy adopted by the cognitive users may be more or less aggressive 
depending on the reliability of their channel sensing. 

Channel identification has a long history in signal processing. The problem becomes especially chal- 
lenging in CR networks, where the estimation of the channel voids, for example, must be very accurate. 
Nevertheless, the estimation itself may be improved by exploiting the availability of a network of nodes 
that could, in principle, cooperate to get better and better estimates of the electromagnetic environment, 
working as a sensor network of cognitive nodes. 
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